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ABSTRACT 

For the purpose of studying the effect of varying quality upon the 
pressure drop accompanying isothermal, two phase, two component flow, 

a horizontal, one-inch flow channel of circular cross section, into which 
air could be injected at various points, was designed and constructed. 

Preliminary tests were made using water and air injected at a single 
position. The resultant data confirmed the findings of Chisholm and 
Laird for turbulent flow in a smooth tube. 

Tests were then conducted in which the quality, the percent mass of 
gas phase of the total mass flowing, was varied by injecting known amounts 
of air at seven injection points along the 42-foot test section. The 
Chisholm-Laird correlation provided pressure drop predictions to within 


25% of the measured values. 
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1, Introduction. 

Accurate prediction of the pressure drop experienced by two phases 
of a substance (or substances) flowing concurrently through a pipe or 
channel is required in many engineering situations. Some common uses of 
two phase flow are: pipeline gathering systems in the gas and oil industry; 
heat exchangers; horizontal tube evaporators; steam generation lines; re- 
actor coolant channels; pipeline contactors in chemical plants; and conden- 
sate lines in refrigeration plants. 

The term, 'two phase flow", used herein, will mean either two compon- 
ent, two phase flow, or single component, two phase flow, according to the 
context. 

As fluids pass through horizontal, closed channels or pipelines, they 
experience a loss of pressure through energy dissipation due to turbulence 
and friction between the walls and the fluid. For single phase flow the 
pressure drop may be predicted with adequate precision using the Darcy- 
Weisbach equation and a friction factor obtained either from the 4oody 
Diagram or, if the pipe is smooth and the flow turbulent, calculated from 
the Blasius formula. Unfortunately, a simple summation of the pressure drops 
due to each of the phases eee separately does not predict the actual loss 
of pressure. 

The increased complexity of computations encountered is due to one or 
more of the following reasons. 

The hydraulic diameter of each of the phases is variable and generally 
undefinable as flow progresses. The cross section available for flow of 
one phase is reduced as the second phase is introduced, The unsteady, in- 


tensely turbulent flow, which is characteristic of high mass flow rates in 





two phase flow, is an efficient energy absorber. Much energy is irrecover- 
ably expended creating ripple, wave, plug, and slug flow. No longer is 
the flow boundary a relatively smooth pipe wall, having become instead a 
rough, mobile, irregular interface between the two phases, In view of such 
problems, considerable effort has been required to obtain good correlations. 
The history of two phase flow investigation extends back to 1939 when 
Boeiter and Kepner /18/* first attempted to obtain a workable correlation. 
In 1942 Dittus and Hildebrand /19/ presented a method which used both mathe- 
matical equations and graphical solutions for the determination of pressure 
drop in oil vapor-oil mixtures flowing through furnace coils. Then in 1944 
came the first generally accepted correlation for the prediction of two 
phase flow pressure drop, that of Martinelli /3/. During succeeding years, 
Martinelli and his associates attacked the problem of pressure drop during 
the forced circulation boiling of water /5/, and produced a revised cor- 
relation in which the liquid fraction (the fraction of the cross-sectional 
area of the pipe occupied by liquid) was integrated into the prediction /4/. 
In 1954 Alves /20/ reported experimental results within 15% of the pre- 
dicted pressure drop, and classified the seven types of flow that were 
observed, White and Huntington /7/ in 1955 proposed a new correlation 
which brought the predicted values to within 12% of the measured values, 
for the very restricted type of flow known as “ripple flow”. In 1955 
Chenoweth and Martin /6/ presented a correlation that came within 20% for 
turbulent flow. Chisholm and Laird /1/ modified Martinelli's flow-rate- 
ratio parameter, introduced the effect of roughness, and obtained a cor- 
relation to within 15% without recourse to a graphical solution. A two 


phase friction factor developed by Bertuzzi, Tek and Poettmann in 1957 /21/ 


Se 
/ / denotes reference in Bibliography. 
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provided the basis for a correlation which produced a standard deviation 
of 20.8%. In more recent years, Hatch and Jacobs /22/ investigated the 

prediction of pressure drop in the two phase flow of hydrogen, and Kordy- 
ban /23/ presented a flow model that provided a logical explanation for 

the applicability of the Chenoweth-Martin correlation /6/. 

A preliminary library search produced no reports which dealt with the 
pressure drop associated with flow of varying quality. Consequently a 
horizontal, one-inch, flow channel of circular cross section, into which air 
could be injected at various points, was designed and constructed. 

Preliminary tests were made using water and air injected at a single 
position. The resultant data confirmed the findings of Chisholm and Laird 
for turbulent flow in a smooth pipe. 

Tests were then conducted in which the quality, the percent mass of 
gas phase of the total mass flowing, was varied by injecting known amounts 
of air at seven injection points along the 42-foot test section. The 
Chisholm-Laird correlation provided pressure drop predictions within 257% 


of the measured values, 





2. Description of Equipment. 

a, Flow Channel Apparatus. The 42-foot test section, Figs. 1 and 2, 
were made of extruded,acrylic plastic tubing of a nominal one-inch inside 
diameter, Seven sections of the tubing, six feet in length, were joined 
in station connector blocks (Fig. 3) using O-ring seals. This construc- 
tion provided very satisfactory sealing with 4 minimum reduction in the 
smoothness of the tube. The actual inside diameter varied from 0.995 to 
1.025 inches both in the diametral plane and along the tube. 

An effective diameter of 1.013 inches was used. By assuming the 
tube to be smooth, using equations (1) through (8) (listed in Section 3-b) 
with the values obtained from the single phase water runs, and treating the 
diameter, D, as the unknown variable, the resultant equation was solved for 
the effective diameter. Thispatancean was then used, as required, in sub- 
sequent calculations. 

Figure 4 displays the results obtained from the single phase water 
runs using this diameter. It is evident from Fig. 5 that the single phase 
air runs suffered no adverse effects from the use of the diameter thus cal- 
culated. Included in the effective diameter were the effects of the eight, 
station connector blocks. Slight discontinuities in the smoothness of the 
tube were introduced by the two tubing ends butting against an insert of the 
same tubing glued inside the block, and by the annular series of holes for 
each of the pressure taps. Distances between stations are shown in Appendix 
1 

b. Auxiliary Equipment. An overall view of the test equipment is pro- 
vided in Fig. 1. A more detailed view of the auxiliary equipment shown in 


Fig. 6. In this photograph the following pieces of equipment are evident: 





Number Item 

1. Water Drain Stop Cock 

2, Rotameter Selector Switch Valve 

3. Rotameter Selector Switch Block 

4, Injection Air Rotameter 

5. Water Rotameter 

6. Differential Mancometer Inlet Pulsation Damper 
7. Station Pressure Manometer Pulsation Damper 
8. Manometer Line Bubble Trap 


9, Water Main Inlet to Bubble Trap (with air bleeder valve 
directly abcve) 


10. Precision Potentiometer 

11. Electric Timer 

12. Timer Control Box 

20. Cathetometer 

21. Manometer Bank for Measuring Rotameter Pressure 

22. Differential Mercury Manometer 

23. Air Pressure Regulators in Series 

24, Differential Water Manometer 
Where these pieces of equipment are indicated on the general apparatus 
schematic, Fig. 7, the numbers correspond. 

The main features of the auxiliary equipment, the selector switches, 
may be seen in Figs. 8 and 9. In these figures most details of their 
construction are apparent. Two details not evident are: 

1. Whenever gluing of the plastic was required, a mixture 
of thin acrylic shaving in ethelyene dichloride was inviscid enough to 


penetrate the thinnest crack or abutment of the pieces being joined. This 





method provided joints of more than adequate strength. 

2. The aluminum alloy (24ST) slide valve, item 2, Fig. 9 
had slots milled along half of its length,3/16 inch steel tubing laid 
therein, and was returned to its cylindrical shape by filling with ep xy 
cement. This eliminated the difficult task of successfully drilling small 
diameter holes for a long distance in a relatively soft material. The holes 
which are evident in the valve, lead directly to the tubing connections at 
the right hand end of the valve. These in turn are connected to the inlet 
and outlet of a rotameter. 

The rotameter selector switch, Fig. 9, allowed insertion of the ro- 
tameter into each station's air supply line successively, only momentarily 
disrupting the continuity of any of the supply lines, and the "steady" state 
of the flow. 

The station pressure selector switch, Fig. 8, was used to transfer the 
manometer pressure leads from one station to the next, and could also be 
used for easily obtaining differential pressures between successive pairs 
of stations. 

The bubble trap, item 8, Fig. 10, kept air bubbles from entering the 
manometer pressure lines. Water was admitted to the top of the bubble trap 
through the valve, item 9. Entrapped air was then bled off the top of the 
trap as it accumulated. Station pressures were measured in all cases by 
mercury or water manometers vented to the atmosphere, All differential 
pressures between stations were measured with a mercury or water differen- 
tial manometer. In the case of differential pressures, all liquid levels 
were measured to within 0.05 millimeters using a Gaertner cathetometer, 

Temperatures were measured using copper-constantan thermocouples refer- 


red to an ice bath. The wiring diagram for this circuit is indicated in 





Fig. 11. Temperature measurements made at the inlet (thermocouple #1) 
and at the outlet (thermocouple #2) indicated that the flow was isother- 
mal to within plus or minus 0.1°F. 

The high pressure air source fluctuated between 80 and 100 psig. 

This fluctuation was reduced to plus or minus 0.1 inches of water by the 
use of two pressure regulators connected in series, 

The water main pressure was found to be sufficiently constant, varying 
in mass flow rate by only 0.007 lbm/sec, when the tests were conducted 
after hours of normal water usage. 

The weigh tank scale was checked with dead weights. 

Rotameter readings were converted to cubic feet per minute using the 
curves supplied by the manufacturer, 

The test section was horizontally aligned, using a K&E Tilting Level, 
to within plus or minus 0.03 inches. A chalk line, snapped between the 
extremities of the test section support, provided a straight reference 
line for longitudinal alignment. See Fig. 2. 

Figs. 12 and 13, with the wiring schematic, Fig. 14, show the method 
employed to measure the time required for a4 set mass of water to flow. 
Movements of the tank switch contact arm caused this mass of water to vary 
from 148 to 154 lbm, and consequently the "set'' mass frequently required 


redetermination. 
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3. Procedure. 

a. General. Except as modified in subsequent sections, each test 
run was conducted in a similar manner. Temperatures were stabilized by 
running both the water and air for an hour before the test. Atmospheric 
temperature and pressure were recorded, The quick closing valve on the 
weigh tank was closed and the weighing proceeded automatically. Water and 
air injection temperatures were recorded, The overall pressure drop, station 
pressures, and air rotameter pressure were then noted. 

b. Single Phase Qualifying Tests. In order to establish the validity 
of measurements made, tests were conducted to determine the friction pres- 
sure drop for single phase water flow and single phase air flow. The re- 
sults of these tests are displayed in Figs. 4 and 5. In all cases in 
this investigation, data reduction was accomplished using a digital com- 
puter, 

The determination of friction factors, Reynolds numbers, and pressure 
drops for the case of a fluid flowing alone was made using the following 


equations: 


(1) LJ = fe Ap Y ; 


where: WA. 


At 


Ap 
VL 


Mass flow rate, lbm/sec, 


Density of liquid, Tenvce: 


i 


iz 
Area of tube cross section, ft , 


Average velocity of flow, ft/sec, the velocity 
of the liquid which would occur if only the 
single phase were flowing; 


rearrangement gives, 
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Using this velocity, the Reynolds number was calculated by: 


Vi Dee 
3) NN. — edbve 
R A, ; 
where: Ne 
D 
fit 


Since the test section was a smooth tube, the Blasius formula provides 


Reynolds number, dimensionless, 


Effective diameter of tube, ft, 


Absolute viscosity of liquid, lbm/ft-sec. 


the applicable friction factor: 


eels : 
(4) = Ne 9 = Friction factor, dimensionless. 
R 


One form of Bernoulli's equation states: 





(5) 2, + Lie + Ye = Z,+ Re 4s Ve aga ? 
Je ‘a Z% Se V4 Za. if 


In the case of horizontal flow with constant area and fluid density, 7 : 
Zz ie aes » and equation (5) reduces to the form: 
(6) i = Me — Af = pA 
Ae jellies 


where: Af 





Pressure drop due to liquid flowing alone in 


the tube, 


Ag = Friction head, ft-lbf£/1bm. 


The Darcy-Weisbach equation states: ; 


(7) dpe bof [LL = Length, ft. 


Equating (6) and (7): 
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where: AR = Pressure drop per unit length, psf/ft. 


c. Two Phase Flow. 
i. Background of the Chisholm and Laird Correlation. 

Use is made in this paper of a two phase flow pressure drop correlation 
developed for smooth and rough tubes by D. Chisholm and A. D. K. Laird 
/1/ in 1956. The reasons for choosing this correlation are twofold; the 
data could be correlated to within 20%, and the portion of the correla- 
tion used, was completely in mathematical terms, which facilitated data 
reduction by a digital computer. 

When correlations which were available to Chisholm and Laird were 
used, the percent error in prediction of pressure drop increased as a 
given set of flow conditions deviated farther from the conditions of flow 
ona smooth tube. These errors were commonly of the magnitude of 30%. 

To produce a better correlation Chisholm and Laird collected data 
from concurrent, horizontal, isothermal flow of water and air through 
smooth and rough tubes of approximately one-inch bore. The ranges of ex- 
perimental data which they correlated, and the types of tubes used, are 


summarized, with those of the present investigator, in Tables 1 and 2 on 


the following pages. 
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TABLE 1 


RANGES OF DATA 





UANT ITY CHISHOLM & LAIRD AUTHOR UNITS 
Mass velocity 9 
water 39 to 600 51 to 236 lbm/sec-ft, 
air O21 to 20 0.589 to 2.10 lbm/sec-ft 
Reynolds number 
water 4500 to 80,000 5790 to 26,000 
air 1000 to 140,000 4090 to 14,580 
Temperature 
water 59 to 74 52.3 to 61.8 degrees F. 
air 64.8 to 68.3 degrees F. 
Viscosity -5 ie -5 ae 
water 77x10_gto 63x10 ~_, 85x10 “te 73x10 _Ubm/ft-sec 
air 1, 2x10" to, azsx10 1.2x10 ~ to 1.21x10 ~1lbm/ft-sec 
Tube Surfaces 
€/), pipe roughness ratio 
Acrylic plastic* 0.000 
Smooth brass** 0.000 
Commercial galvanized tube** 0.0025 


Brass tube, concrete internal surface** 0.013 


Brass tube, internal thread** 0.028 (measured) 
0.037 (apparent) 


Brass tube, non-uniformly distributed sand**0.045 
Brass tube, uniformly distributed sand** 0,068 
General Remarks 

Flow in both phases was considered turbulent throughout the ranges 
considered, Mean arithmetic values of the water temperature at the inlet 
and the outlet were used when determining water properties. Air properties 
were determined with no appreciable error by assuming that the air was dry 
and the temperature quickly became the same as the mean water temperature 
upon mixing. 


* Author's tube 
** Chisholm & Laird's tube 11 





TYPE OF FLOW 


TABLE 2 
SPECIFIC RANGES OF MASS FLOW RATE 
FOR EACH TYPE OF FLOW 


STUDIED IN THIS INVESTIGATION 


WATER, lbm/sec AIR, lbm/sec 


Single phase water 0.508 to 1,32 = ee 


Single phase air 


~~ - 0.00330 to 0.0118 


Two phase single injection 0.286 to 1.27 0.00348 to 0.00950 
Two phase multiple 0.694 to 1.30 0.000426 to 0.00142 
injection (per injection point) 
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Through introduction of the effects of roughness into the Martinelli 
correlation Chisholm and Laird arrived at the following correlation of 


their data. They define a flow-rate-ratio parameter, x » by equation 
675 a I25 Ay 


x= @ A) & 


Liquid-mass velocity based on tube cross section, temic ce 


where; Ge 
Ge = Gas-mass velocity based on tube cross section, meee tee 
Mi= Absolute viscosity of liquid, 1bm/sec-ft, 
As = Absolute viscosity of gas, lbm/sec-ft, 


St 


Je = Density of gas, Meee 


Density of liquid, lbm/fe°, 


| 


=D 


The pressure-drop parameter, Aftp 9 is then correlated by equation 


A 


AFre 
AR 





(10) 


where: AfPp= Total (friction) pressure drop over an increment of length 


Are = Friction pressure drop for liquid flowing aloné in tube, psf/ft, 


C=> War 


M= 


= Constants for a particular liquid flow rate and tube surface, 


»>|>- 


rp= Reynolds number where the liquid flows alone, 
» = Pipe friction factor for rough tube, 
dg = Pipe friction factor for smooth tube. 
Equation (10) may now be expressed as 


B7Sim (25m 


(11) Sie = | +c(@ a 
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The liquid fraction, RRL » is correlated with the pressure-drop para- 


meter, APre » as follows: 


AR 


(12) a = ~~ (for the smooth tube), 





where: KR. «= Liquid fraction, the area occupied by the liquid 
divided by the cross-sectional area of the tube. 


A correlation of the liquid fraction, Ree » with the flow-rate-ratio para- 


meter, xX » for the smooth tube is; 


(13) 08. = | + 2! — a (for the smooth tube). 


Rie a x x 2 


ii. Single Injection Tests. After the single phase tests were 


completed, a series of runs was made in which water and air were injected 
upstream of the test section. The correlated pressure drop was obtained 
in the following manner: 

1s X was determined using equation (9). 


2. Rt was obtained from equation (13). 


ci AF cout then be calculated from equation (12). 

The results of these single injection tests, together with the results 
of the multiple injection tests (described in the next section), are pre- 
sented graphically in Figs. 15, 16, and 17, where the coordinates used are 
those found to be the most useful by Chisholm and Laird. 


iii. Multiple Injection Tests. Having thus shown that the ap- 


paratus is capable of producing data comparable to that presented by 
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Chisholm and Laird, test runs were conducted in which air was injected 
in known amounts at seven points along the test section. 

The purpose of this multiple injection of air was to provide portions 
of the test section with two phase flow whose liquid fraction was able to 
be determined. In addition, it was desired that the liquid fraction vary 
from injection point to injection point. The design of the test section 
allowed injection of air at each station connector block. Since the amount 
of air injected at each station was measured, the liquid fraction was calcu- 
able, and constant for each portion of the test section between injection 
points. The correlation process, set forth previously in section 3-c-ii, 
was then used and the pressure drop found for each constant-liquid-fraction 
section. The pressure drops were then summed over the whole test section, 
and the overall pressure drop obtained. The results of these multiple in- 


jection tests are displayed in Figs. 15, 16, and 17. 


se) 





4, Results and Conclusions. 

A horizontal channel suitable for multiple air injection into the 
primary flow of water was designed, constructed, and tested. 

The results of single phase tests, Figs. 4 and 5, indicated that 
the test section behaved as a smooth tube. Consequently, the friction 
factor computed by the Blasius formula was the only one employed during 
subsequent reduction of the two phase runs. 

These Figures show that the pressure drop accompanying two phase flow 
of varying quality may be predicted to the same degree of accuracy as 
constant quality two phase flow, by a suitable selection of the "constant 
quality" sections. A logical next step would be the prediction of pressure 
drop in a two phase, varying quality, isothermal situation, such as a 
Saturated steam generating tube, by applying the correlation to sections 
of the tube small enough to be of "constant quality", from one end to the 
other. This could be done on a digital computer and the results compared 
with experimental data. A further extension of this line of investigation 


might well start with mounting the variable quality test section vertically. 
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APPENDIX I 


LOCATION OF PRESSURE TAPS ALONG THE TEST SECTION 


Station Distance to next Distance to station indicated 
Number station,in inches from station one,in inches 

1 71.7 0.0 

2 WVe72 vod 

) 71.75 143.42 

4 71.8 215.17 

5 71.75 286.97 

6 VES 358.72 

7 71.7 430.52 

8 0.0 3025 22 


Distance of station 1 from inlet --37.22 inches 
Distance of station 8 from outlet --35.53 inches 
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APPENDIX IT 


DRAWINGS, PHOTOGRAPHS, AND GRAPHS 
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Fig. 1 General View of Equipment 








2 Overall View of Test Section 
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Fig. 4 Friction Factor vs. Reynolds Number 
Single Phase Water Flow. 
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